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Abstract. 

We derive an expression for the entropy of a dark matter halo described by a 
Navarro-Frenk- White model with a core. The comparison of this entropy with the one 
of dark matter at the freeze-out era allows us to constraint the parameter space in 
mSUGRA models. Moreover, joining these constraints with the ones obtained from 
the usual abundance criteria and demanding both criteria to be consistent with the 
2a bounds for the abundance of dark matter: 0.112 < ^lr>Mh 2 < 0.122, we are able to 
clearly discriminate validity regions among the values of tan/3, one of the parameters of 
the mSUGRA model. We found that for the explored regions of the parameter space, 
small values of tan/3 are not favored; only for tan/3 ~ 50 are both criteria significantly 
consistent. In the region where both criteria are consistent we also found a lower bound 
for the neutralino mass, m x > 141 GeV. 
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1. Introduction 



The existence of components of the total energy- density of the Universe whose nature is 
still not understood, constitutes one of the biggest unsolved questions in physics today. 
It is a great challenge to have a clear picture of the nature of dark matter (DM), and dark 
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energy, whose existence is more and more undisputed with the observational evidence 
accumulated over the last decades. 

The observational constraints on the present density of DM, Qdm, come from 
several outstanding observations such as the Cosmic Microwave Background radiation 
(CMBR) P], galaxy clustering, supernovae and Lyman a forest. One of the most recent 
works which combines all these data leads to: 0.112 < Q^Mh 2 < 0.122 [2], with h the 
dimensionless value of the present day Hubble constant (h ~ 0.7). These bounds imply 
that dark matter constitutes 22.8 — 24.8 % of the total density of the Universe. 

One of the most accepted candidates to be the major component of dark matter 
is the neutralino as the Lightest Supersymmetric Particle (LSP). Supersymmetric 
theories with R-parity conservations predict this particle (for an excellent introduction 
to Supersymmetry see [3]). This type of models have several parameters, some of which 
can be constrained by combining different methods with observational data of the actual 
density of DM. In particular for mSUGRA models, this has been done using a standard 
approach [U [5] , based on solving the Boltzmann equation by considering that after the 
"freeze-out" era, neutralinos cease to annihilate keeping its number constant. In such 
approach, the relic density of neutralinos is approximately given by (this is strictly true 
only when (av) is independent of the energy): 

n x oc l/(av), (1) 

where (av) is the thermally averaged cross section times the Moller velocity of the 
annihilating pair. 

Within the mSUGRA model, five parameters are needed to specify the 
supersymmetric spectrum of particles and the final relic density of the LSP. These 
parameters are: m , the unified mass for scalars, m^, the unified mass for gauginos, 
Aq, the unified trilinear coupling, tan/?, the ratio of the vacuum expectation values of 
the neutral components of the supersymmetric Higgs bosons and the sign of fi, where \i 
is the Higgsino mass parameter. In this work we use the numerical code micrOMEGAs 
[6] to compute the LSP's relic density for different values of these parameters. This 
scheme will be called hereafter the "abundance criterion" (AC), which combined with 
the hypothesis Qdm = ^x-> gi yes an effective way to constraint the values of those five 
parameters. 

Another estimate for Q x can be obtained using a different approach. Just before 
the "freeze-out" epoch, we can consider neutralinos as forming a Maxwell-Boltzmann 
(MB) gas in thermal equilibrium with the rest of the components of the Universe. In 
the present time, such a gas is almost collisionless and either constitutes galactic halos 
and larger structures or it is in the process of their formation. In this context, we can 
conceive two equilibrium states for the neutralino gas, the decoupling (or "freeze-out" ) 
epoch and its present state as a virialized system. Computing the entropy per particle 
for each one of these states, we can use an "entropy consistency" criterion (EC), based 
on the comparison of theoretical and empirical estimates for this entropy, to obtain an 
alternative expression for the relic density of neutralinos Q x , functionally dependent of 
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the mSUGRA parameters. This idea was originally introduced by some of us in [7|. 

Our objective is then to develop the recently introduced EC, as a complementary 
method to the well known AC, in order to obtain further constraints on the parameters 
of the mSUGRA model by demanding consistency of both criteria with each other and 
with the observational constraints on Qdm (preliminary results of this analysis were 
already presented in [8]). 

The paper is organized as follows, in section II we present a brief description of the 
method to compute fl x using the AC. In section III a derivation of the EC is explained. 
In section IV we present the halo model together with a detailed description of the 
method we followed to obtain one of the key empirical parameters in the EC. In section 
V we present the results obtained by comparing both criteria with the observational 
constraints on Qdm in the context of the mSUGRA model. In the last section a summary 
and the conclusions of our work are presented. 

2. Abundance criterion 

The relic abundance of neutralinos is defined as fl x = p x j p cr it, where p x = m x n x is the 
relic's mass density (n x is the number density, m x the neutralino mass) and p cr i t is the 
critical density of the Universe (see [9] for a review on the following method to compute 
the relic density). The time evolution of n x is given by the Boltzmann equation: 

^ = -3Hn x -(av)(n x -(n x «n (2) 

where H is the Hubble expansion rate and n e x q is the number density of neutralinos in 
thermal equilibrium. In the early Universe, neutralinos and the rest of species were 
in thermal equilibrium, that is n x = n e x q . As the Universe expanded, their typical 
interaction rate started to diminish and the process of annihilation froze out. Since 
then, the global comoving number density of neutralinos has remained nearly constant. 

There are several ways to solve Eq. §2§. The most common approach is based on 
the "freeze-out" approximation (see for example |J0]). If we consider that neutralinos 
are the major component of the dark matter density, Q x ~ Qdm, then the freeze- 
out approximation gives a functional dependence between Qdm and the mSUGRA 
parameters. In other words, it gives a constraint equation for these parameters. 

In practice, however, we obtain a numerical solution to the Boltzmann equation, 
Eq. ([2]), using the public code micrOMEGAs 1.3.6 [6], this gives us better precision 
over the usual freeze-out approximation. The code computes the relic density of the 
LSP in the Minimal Supersymmetric Standard Model (MSSM), taking into account all 
annihilation and coannihilation processes as well as loop-corrected masses and mixings to 
calculate (av) exactly. Further on, we take the mSUGRA model and its five parameters 
(m , m 1 / 2 , A , tan/3 and the sign of /i) as input for micrOMEGAs and use Suspect in 
its version 2.34 which comes as an interface to micrOMEGAs, to calculate the 
supersymmetric mass spectrum of particles. 
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Using this criterion, we can obtain the allowed regions in the parameter space of 
the mSUGRA model which are consistent with observational constraints on the actual 
dark matter density. 

3. Entropy consistency criterion 

This criterion was originally introduced in [7J. We briefly describe it in the following. 

The main idea of the method is to compare the neutralino gas in two stages of its 
evolution as the Universe evolved, namely, the freeze-out era and the present epoch. 
These initial and final states are taken as equilibrium states. In this approach we 
disregard the complex phenomena that took place during the formation and evolution 
of structure in the Universe, and the analysis is focused only on those two stages. In the 
freeze-out era, the assumption is to treat the neutralino gas as a non-relativistic ideal 
gas of WIMPs (weakly interacting massive particles) described by Maxwell-Boltzmann 
statistics. This assumption is valid mainly because in this era, the dominant interactions 
in the gas are still short-ranged and neutralinos are already non-relativistic by then 



We would like to apply the MB description for present day structures as well, in 
particular to galactic halos, which we consider as final states of the primordial neutralino 
gas. However, the valid assumptions made for the freeze-out era are no longer valid for 
such final states, because in galactic halos, neutralinos are mainly subject to long- 
range interactions, related to non-extensive forms of energy and entropy. Following 
the guidelines given in [7], a more general approach can be applied instead, it uses the 
microcanonical ensemble in the "mean field" approximation. This approach is actually 
valid at both, the initial and final stages. The state of the system in each state is defined 
by its entropy per particle s, its density n, and its temperature T; the latter is actually 
given by the auxiliary variable x — This set of variables is then (s^, x^, nt) for the 
neutralino gas at the freeze-out era and (s (h) , x (h) , n (h) ) for the dark matter halo today. 
Therefore, the change in entropy per particle between these two states is given by: 



The region today in which all these considerations apply is the center of galactic halos. 
Thus, in what follows, we consider current halo macroscopic variables as evaluated in 
this region: s^, xf \ n^ } . 

We can then write s^ using Eq. (j3J) as an equation depending only on present day 
values of certain quantities, such as the ratio of the total density of the Universe to p cr it, 
Qq, and the value of the freeze-out temperature of the neutralino gas. This was done in 
detail in [7] and the result can be written as: 



{m x > 100 GeV). 




(3) 




th 



- + Xf + In 
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where g. t is related to the degrees of freedom of the system, Xq mb = m x /T^ MB = 
4.29 x 10 12 m x /GeV, with T CMB = 2.7K (see \10\ for details). 

It is important to mention that there is an alternative formalism to study 
autogravitating systems such as dark matter halos. It is based on a non-extensive 
description of energy and entropy developed by Tsallis p2] . Such formalism was applied 
to dark matter halos and compared with the predominant ACDM paradigm obtaining 
suggestive results [TSj, IT4"]. However, for the purposes of this work we will continue with 
the traditional formalism. 

If the neutralino gas in present halo structures would strictly satisfy MB statistics, 
the entropy per particle would follow from the well known Sackur-Tetrode entropy 
formula [15] . Such a MB gas in equilibrium would be equivalent to an isothermal halo 
if we make the connection [16]: <r 2 h) = k B T {h) /m, where <r ( 2 h) is the velocity dispersion (a 
constant for isothermal halos). However, as mentioned before, the assumption of MB 
statistics does not apply for self-gravitating collisionless systems. Hence, an exactly 
isothermal halo is not a realistic model, not only because of these theoretical arguments, 
but also because its total mass diverges and its distribution function allows for infinite 
particle velocities (theoretically accessible in the velocity range of the MB distribution). 

More realistic halo models follow from "energy truncated" distribution functions 
[16]-[20j that assume a maximal "cut off" velocity (an escape velocity). Therefore, 
we can provide a convenient empirical estimate of the halo entropy, sj, h) , following 
the microcanonical definition of entropy in terms of the allowed phase space volume, 
restricted to the range of velocities accessible to the central particles, that is up to a 
maximal escape velocity v e (0). If we reasonably assume that: 

v 2 e (0) = 2|$(0)| ~ aaf h) (0), (5) 

where $(r) is the Newtonian gravitational potential, and a is a proportionality constant, 
then we can write: 

/ m \ 4 / a \ 3 ^ 2 GeV/cm 3 



sl h) | ~ln 



c 



4 

m v% 



(h) 



89.17+ln 



VGeW 



J'C / fJc 



(6) 



where we used x^ = c 2 /<r ( 2 h) (0) with c the speed of light. 

Equating the theoretical and empirical estimates for the entropy per particle (Eqs. 

and ([6])) we finally obtain a relation for the relic abundance of neutralinos (this 
formula is a small modification of the one presented in [7]): 

{x { a) 3 / 2 m x g, {x% MB 



\n(tt x h 2 ) = 10.853 -x { + In 



(7) 



9*f (*f) 

In this way we have obtained another constraint equation for the parameters of the 
mSUGRA model (recall that VL DM = Q x ). 

In order to perform an analysis of the mSUGRA parameters using this new entropy 
criterion, we modified micrOMEGAs to compute the value of £f for any region of the 
parameter space, and then compute Q x using Eq. ([7]). We can then find the subspace 
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of values of the mSUGRA parameters consistent with the observed constraints on the 
dark matter density. 

The extra parameter a in Eq. (CO), which is the proportionality factor between the 
escape and dispersion velocities at the halo center has a major role in our work. It is 
very important for the EC that we are able to give values to a in an independent way. 
In the following section, we show in detail how this parameter can be evaluated for a 
particular dark matter halo model. 

The scheme we have developed so far is only for dark matter particles, namely 
neutralinos. We have not included at any moment the baryons which constitute the 
galaxies that we observe today. As we discuss further on, this condition makes our 
method strictly valid only for pure dark matter structures and gives a reasonable 
approximation when the dynamical effects of baryons can be neglected. The main 
effects coming from the addition of baryons that should be included in our model are: 
i) neutralino-nucleon elastic scattering interactions and ii) gravitational effects from the 
baryonic component onto the dark matter halo. Although both of these effects are 
more significant in the central region of the halos, precisely where our method should 
be applied, we believe that the latter can be avoided if we focus our analysis on systems 
where the baryon-to-dark mass content is significantly low, and the former can be safely 
neglected in all cases since the interaction rate between neutralinos and baryons is highly 
suppressed by the low neutralino-proton scattering cross section predicted by theoretical 
models 10~ 44 — 10~ 45 cm 2 (see for example [38] , see also [39] for recent experimental 
upper limits). The formation of a baryonic component in the galactic (halo+baryons) 
system, a disc for example, modifies the original mass distribution of dark matter, the 
effect is a net contraction of the halo towards the center, see for example [44] for a 
classical description on how to include this effect. In what follows we will not consider 
it though, since we will focus on dark matter dominated systems where the effect is less 
significant. 

4. About the parameter a 

The formula for the entropy per particle in the center of halos, Eq. ([6]), is far from being 
the final description from a statistical-mechanics point of view. The question about 
the actual description governing "dark matter fluids" remains open. Nevertheless, it 
is the assumption we have made about the neutralino gas, taking it nearly as an ideal 
gas, that make Eq. ([6]) consistent. In fact, other works in the past have considered the 
approximation of dark matter as an ideal gas as the approach to follow when it comes 
to analyzing the entropy (see for example [40], |41|). Our ignorance on the correct 
statistical-mechanics treatment for systems formed by dark matter is reflected in the 
appearance of the parameter a in Eq. ([6]). It is of key importance to find appropriate 
bounds for its value. In the following, we present an specific model used to obtain 
estimates for such bounds. We have chosen this model trying to balance simplicity with 
an approximate general description of dark matter halos. 
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The model consists of an spherical dark matter halo with a constant density core in 
the center, beyond which the dark matter density profile follows the well known Navarro- 
Frenk- White (NFW) profile [12], we set a cut-off radius to define the halo boundary: 

p c if y < y c 



i>i>D= { ^ if y c <y<y v 

if y > y v 



(8) 



where p c is the constant central density of the core, y = r/r s , y c = r c /r s , y v = r v /r s , 
r s is a scale radius, r c is the core radius and r v is the virial radius, which defines the 
limit for the virialized structure; po, 5q and r s are the usual parameters of the the NFW 
profile. The choice of an NFW profile for the external part is clearly motivated by the 
success of the profile as a universal fit to dark matter halos in numerical simulations. 
The central core is motivated by the apparent need of it in dark matter dominated 
systems, see further below, and also by the phenomenon of dark matter annihilation 
which could produce it naturally (see for example [43J ) . 

For a model without core, all these parameters can be given by a series of well- 
established formulas (see for example [Til H2l PH I H5]): 

Ac 3 

6 ° = 3[ln(l + c)-c/(l + c)]' (9) 

Po = p c <Ah 2 = 277.8 h 2 ^, (10) 

kpc 

where c = r s /r v is the concentration parameter. We choose a ACDM model with 
Q = 1, Q\ = 0.7, which seems to be in accordance with different sources of observations 
[JJ; and A ~ 100 061 PJ7]. 

With Eqs. ([9]) and (flOl) and the definition of the parameter c, the NFW density 
profile is defined by two parameters only, a "size" parameter (r s or r v ) and the 
concentration parameter c. In fact, both can be written in terms of the total mass 
of the halo, M v : 

1 3M " V' 3 „„ 



co- 160 . (12) 

Eq. (Tl2|) is a fit extracted from the numerical study [H], which is based on the 
Millennium simulation data pJ2], the sample of relaxed halos presented in that work 
was chosen as more appropriate for our purposes (the halos in this sample are closer to 
virial equilibrium). The results were transformed from the original value A = 200 to 
A = 100, to do so, we used the code provided by the authors. Other studies in the past 
have found different fitting formulas for the mass-concentration relation, we believe our 
results are not highly sensitive to a specific formula (we actually tested this with the 
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formula used in[50j). Eq. (1121) however, just gives a central value for the concentration 
(this is why we use the term cq to represent it), actually, the relation between c and 
M v in numerical studies has a scatter, ai ogc , which is slightly dependent on M v itself. 
For the remainder of this work, we take however the same value for all halos, given by 
the mean of o~i ogc for the sample of relaxed halos presented in [18] (see table 1 of their 
paper): a logc « 0.095. 

Using Eqs. fl9l fT2l together with the NFW density profile, without core, we have a 
model that can in principle be used to describe dark matter structures with only one free 
parameter: M v . A particular case can always deviate from this average model mostly 
because of the scatter in the concentration parameter c. This situation can be alleviated 
if one takes into account the expected range of values for c given by oi ogc « 0.095 as 
described above. 

Such description is, however, not valid for all dark matter halos. There is specially 
a controversy, still not resolved, of whether or not all galaxies are consistent with the 
cuspy behavior of the NFW profile in the center. For dwarf and low surface brightness 
(LSB) galaxies there seems to be a disagreement |51]-|56j. On the other hand, for 
large structures like galaxy clusters, the NFW description is strongly supported by 
observations [57]- [59]. Recent numerical simulations are able to resolve regions closer 
to the center of dark matter halos and seem to favor even cuspier density profiles. For 
example in [60J, p ~ r~ 7 with 7 ~ 1.1, 1.2 and 1.35 for cluster-, galaxy- and dwarf- 
size halos respectively, another analysis ([EE]) finds 7 ~ 1.2 for a cluster-size halo. 
Both analysis agree that a model with a core underestimates the simulated dark matter 
density within the resolution limits of each work, this limit is roughly 1 percent of the 
virial radius. The estimates for the parameter a that we obtain further on are still 
consistent with these numerical results since the values for the core radius that we find 
for our model are slightly inside this resolution limit. 

For the purposes of our analysis, and given the hypothesis we have made, we are 
interested in structures that are dominated by dark matter, that is, that the baryonic 
or visible component is less representative and has no important global effects. In fact, 
in order to empirically estimate the entropy of dark matter halos, we have neglected 
the effects of the luminous galaxies within them, thus, our model should be in principle 
formally valid only for structures that are made exclusively of dark matter. Dark matter 
dominated systems, like dwarf and LSB galaxies, whose dynamics is close to systems 
which only have dark matter, will be taken as a reference to put bounds on the parameter 
a. Taking into account the information we have on such systems, the choice of a halo 
model with the same properties of an NFW profile in the external region, but with a 
core in the central region, seems to be reasonable. 

We consider Eqs. (l9"]fT2l) to be valid as well for this model. Such assumption is 
strictly valid for Eqs. (l9lfT0l) because they refer only to a particular way to define the 
normalization of the profile, they can be used as long as one takes into account its 
connection to the central density: Sopo/(y c (l + y c ) 2 ) = Pc- The same reasoning is valid 
for Eq. (TTTj) . which is independent of the particular model for the density profile, it is 
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a definition for virialized spherical systems depending only on the cosmological model. 
Eq. (j!2p is valid only as an approximation, the virial mass appearing on it refers strictly 
to the virial mass computed for an NFW model without core, the introduction of the 
core reduces the mass in a percentage that depends on how large is r c . However, for 
low values of r c compared to r v , the difference caused by the core is not relevant and 
the formula can be taken as a good approximation. In any case, Eq. ffT2l can be always 
considered as a fair way to give a value to one of the free parameters (c or r s ) in the 
model. 

With all the previous considerations, we can continue and calculate analytically the 
parameter a. First, the mass profile for our model follows from Eq. ([8]): 



M(y) 



^p c r 3 s y 3 if y < y c 



Wo 2 {im + ln(^)-^) Hy c <y<y v 

(13) 



where Vq = 4irGr 2 5 p and G = 4.297 x lO _6 (km/s) 2 kpc/M is the gravitational 
constant in appropriate units (here we have followed closely the work [62]). Using 
the mass profile, we calculate the gravitational potential at r = 0, $(0): d§(r)/dr = 
GM(r)/r 2 ; we use also the fact that $(r„) = —GM(r v )/r v (p = for r > r v ), we obtain: 

™-«{Tk-w&Y (l4) 

In order to calculate the velocity dispersion in the center, cr(0), we use the Jeans equation 
for spherical systems [161 EH] : 

l d(pv 2 ) | 2 Bv? _ d<$> 
p dr r dr } 

where v 2 is the mean radial square velocity and B = 1 — v^/v 2 , with v 2 e the mean 
square velocity in the 9 direction. B is a measure of the anisotropy in the system. It is 
worth saying that several studies have shown that dark matter structures are actually 
anisotropic (see for instance [63| 164"]); however, in the center of these structures, the 
studies usually give B ~ 0. We stay within the isotropic model for simplicity considering 
that it will be a good approximation in the central regions to real structures with 
a density core. Such assumption implies also that v 2 = a 2 , where a is the velocity 
dispersion, related to the temperature of the neutralino gas. We solve then the Jeans 
equation for the two regions in which our profile is divided (y < y c and y > y c ) and 
match both solutions (recalling that p{r v ) = a(r v ) = 0). Thus, the velocity dispersion 
in the center is: 

^(0) = VO'VCCW {6(jL ^ (l + _ Lh{l + yc)) + 



lnfi±^ 
-lnf^ 



qlr, / l±Vv\ _L 2y e (% e +3) , j/g-5i/g-3y„+l 

0111 1 i+uc ) + (i+j/c) 2 + vW+yv) 



61., (-^) • 1 • 2 



+ 1 



lh $fi ( -23y, 3 - 41^ 2 - 9y v + 3 + ^ (-12^ - 18y 2 



yv(i+yv) 
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4-1. Analysis on LSB galaxies 

Using Eqs. (1141) and (I16p we have an analytical expression for the parameter a. In order 
to give an empirical estimate for the range of values that a can take, it is necessary to 
compute Vq, y c and y v , which is equivalent to give values for r v , r s and r c . To do so 
we use observational data extracted from the recent sample of LSB galaxies presented 
by [65]. Although LSB galaxies are dark matter dominated systems, baryons have non- 
negligible effects in the total mass distribution of the galactic system, specially in the 
central region. A more realistic approach to these systems consists then in taking into 
account the dynamical effects of baryons as an extra gravitational component, which 
changes the original dark matter distribution, pulling it inwards towards the center 
during the process of disk formation (see for example [44] for a description of a classical 
model that incorporates this effect). Such analysis is more complicated and we consider 
that it would not increase significantly the accuracy of our model; the main reason is 
that the formula for the entropy of halos (Eq. (jSJ)) was postulated under the assumption 
of a single dark matter fluid, the addition of a baryonic component could invalidate it, 
the mass distribution model should then be changed only if the entropy formula changes 
accordingly. We believe then that taking the so called "minimum-disk" approximation 
(or "zero-disk"), that is ignoring the contribution of stars and gas, is consistent with the 
hypothesis we have used so far. As we emphasize later on, the values of the parameter 
a found using this approximation on the sample of LSB galaxies, are roughly consistent 
with a particular analysis made for a more clear prototype of a dark matter dominated 
galaxy, where the contribution of baryons can be safely ignored. 

The sample presented in [65] consists of 11 LSB galaxies. We use their table 1 to 
fit our halo model so that it has the same central density in the core, p c (computed by 
the authors using a pseudoisothermal halo model), and finding the best fit to the whole 
rotation curve of each galaxy in the sample. Following the model we have proposed, 
we find a unique fit to each galaxy for each value of the concentration parameter. 
To complete the analysis, we take the expected dispersion on the concentration value 
{o~iogc ~ 0.095) and make fits to the sample of galaxies using the extreme values of 
concentration given by this interval. 

As a result of this analysis, we find a dynamical range in virial masses for the halo fits 
of: 4 x 10 10 M o < M v < 3 x 1O 13 M , or equivalently: 89 kpc < r v < 810 kpc; the range 
of values for the concentration parameter goes from 6.6 to 19.4, and finally, for the core 
radius, we find 0.3 kpc < r c < 4.3 kpc. If we divide the fitted halos in three virial mass 
bins then the following intervals for a are found: 15.5 < a < 32.9, 12.7 < a < 106.4 and 
17.8 < a < 94.9 for the low-, mid- and high-mass halos respectively. Since we would 
like to have an interval for a which is scale independent, we take the final bounds for a 




(16) 



where the dilogarithmic function is defined as: Li 2 (x) = ^4dt. 
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to be: 

17.8 < a < 32.9 (17) 

Of the 11 galaxies in the sample, 7 can be fitted very well, 3 reasonably well and only 1 
poorly with values within this range. In a preliminary analysis to obtain the final interval 
for a, another sample of galaxies was also analyzed, that sample consisted of 6 dwarf, 9 
LSB and 2 low-luminosity galaxies compiled in [66]. For that sample we used a slightly 
different method to obtain a, the two main differences are: i) a different formula for the 
central value of the concentration and its scatter (see [50]) and ii) instead of fitting the 
whole rotation curve, only the value of the maximum rotational velocity was chosen to 
be the same; the interval for a was found to be [16.4,27.8], which is similar to the one 
in Eq. ([171) . Such similarity is found also with the work that motivated the present one 
[7], the authors used a different halo model and followed a completely different path to 
obtain the range of values of a, they found: 11.2 < a < 24.8. As can be seen, despite 
the difference in the halo models and in the methods to obtain a, all these analysis end 
up with similar intervals for its value. 

We consider that the method we have described so far is well suited to give an 
empirical interval for a, which is at the end one of the key parameters in our work. It 
is however important to say that a is in principle a model dependent quantity whose 
value could change by assuming other halo models, or by incorporating other effects, for 
example the addition of the baryonic component or the inclusion of anisotropy. As have 
been said before, the formula for the entropy of halos (Eq. ([6])), which is the basis for 
the final constraints we have found on the mSUGRA parameter space (see next section), 
is strictly valid only for dark matter halos without baryons. So a better alternative to 
get closer to this condition is to analyze dwarf spheroidal (dSph) galaxies, which are 
dark matter dominated systems with very high mass-to-light ratios (higher than a 100 
Mq/L & ). To extend our analysis further and be more certain on the numerical values 
of a, we analyze a specific case, the Draco dSph galaxy (hereafter, Draco). 

4-2. A particular case analysis for dSph galaxies: Draco 

Draco is a good generic example of dSph galaxies since it has been the subject of 
extensive observational studies (see [67] for a recent one). The observational quantity 
measured for these type of systems is the line-of-sight velocity dispersion of its stars, 
ai os , as a function of their projected radius. We use the observational values reported 
in [68] for ai os (extracted from the upper panel of their fig. 5). They took the original 
observational stars sample of [67] and removed unbound stars following a rigorous 
method, see [68] for details. The observational data appears on figure [Q 

To compute c% s for our halo model, we follow the formulas described in [69] and 
[70] . Of most importance is the use of the integral formula for the line-of-sight velocity 
dispersion presented in the appendix of [70] (isotropic case, B = 0): 
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Figure 1. The line-of-sight star velocity dispersion as a function of the projected 
radius for Draco. Observations were taken from [68] . The solid and dotted lines are 
fits from our halo model, see text for details. 

where £(x) = x _p exp(— x l / m ) is the dimensionless 3D luminosity density of the stellar 
component (which can be obtained by deprojecting the surface brightness profile, 
modeled in \69\ by a Sersic law), p = 1 - 0.6097/m + 0.05463/m 2 ; M(x) = M DM (x) + 
f*M*(x) is the dimensionless mass profile consisting of dark matter, M DM , and stellar, 
M*, components; a s and m are the free parameters of the Sersic profile, and finally 
x = r/a s , X = R/a s . Eq. ffT8]) is generically valid for any halo model with spherical 
symmetry. The stellar component has 3 parameters: a s , m and /*, the first two are 
constrained by the brightness profile of Draco [71] : a s = 7.3 arcmin (~ 0.17 kpc for a 
distance to Draco of 80 kpc) and m = 0.83. And /* = M*/M v , where M* = 6.6 x 1O 5 M 

Using Eq. fflBl together with our halo model, we can find fits to the observational 
data using 2 free parameters: M v and p c , and using the central value for the 
concentration parameter (Eq. (fT2"]l ) altogether with its expected dispersion o\ ogc ~ 0.095. 
In Fig. (00) we show different fits to the observational data, the black (thicker) and 
red (thinner) lines are for the upper and lower values of the concentration parameter 
respectively. We found that for the lower values of the concentration parameter, red 
(thinner) lines, the best fit to the data is a model which converges towards a halo with 
no core (red (thinner) dashed line). Increasing the core radius reduces the goodness of 
the fit, we measured this evaluating x 2 /df, where olf is the number of degrees of freedom: 
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number of data points - number of parameters. On the contrary, for higher values of 
the concentration parameter, black (thicker) lines, the best fit is found for a model 
with core (black (thicker) solid line), the goodness of the fit is statistically reduced with 
decreasing core radii, the convergence to a halo with no core is shown as a dashed black 
(thicker) line. These results imply that our model with high values of the concentration 
parameter is able to give an adequate description of the actual observations for Draco. 

Our whole model is built on the assumption of a model with core, this goes back to 
Eqs. P| and (EJ), so we should clearly concentrate on the fits corresponding to a model 
with a core. In the search for the best fit to <j/ os , this can be done unambiguously only for 
certain values of the concentration (those on the high end of the expected concentration 
interval). For instance, the red (thinner) solid line is as good fit to the observational 
data as the black (thicker) solid line, however, we can not extract the value of a without 
ambiguity for the latter since reducing the core radius increases the goodness of the 
fit, recall that convergence is reached for r c — > kpc. The core model which we find 
to be a good fit to the observational data (black (thicker) solid line) has the following 
parameters: M v = 1.5 x 10 9 M o , p c = 3.1 x 10 8 M o kpc 3 , r c = O.llkpc, r v = 30kpc, 
c = 26.6 and a = 28.7. 

We believe that this value of a for Draco is reasonable and certainly consistent 
with the interval obtained before (Eq. (fl7|) ). As a complementary check, we have used 
another set of data on the velocity dispersion of Draco ([72]) to fit our model. Although 
larger values for the virial mass and core radius are needed to achieve a good fit, we 
found the value of a to be consistent with the interval of values given in (Eq. ( flTj) . 
Nevertheless it should be stressed that better observational constraints are needed to 
put firmer constraints on the value of a, which is particularly sensitive to the value of 
p c (r c ). The value of a is of course also sensitive to the chain of assumptions needed to 
arrive to its final estimate, concentration values have an important role as can be seen on 
figure [TJ The assumption of isotropic velocity dispersion, B=0 in Eq. ( fl5i) . has also an 
influence in the values found for the model parameters in the fit to Draco. The addition 
of anisotropy to our model would alter the values of the best fit parameters found for 
Draco, and therefore the value of a we just reported. The impact of anisotropy on the 
value of a can only be properly quantified by improving our model with the removal of 
the restriction B=0, and redoing the analysis we have made so far. Such task is left as a 
possible future work once further more precise constraints on the dynamical properties of 
dark matter dominated systems like Draco are available. We believe that the assumption 
of isotropy is sufficient for the purposes of this work to calculate approximately the value 
of a. 

Overall, taking into account all the assumptions we have made to fit our model with 
observational data, we consider that the interval of values obtained for the "ignorance 
parameter" a, do give a reasonably description to the several galactic samples for which 
we have tested our model, thus making us feel confident on the overall entropy criteria 
presented in this work, at least as a first order indirect analysis of real dark matter 
structures. 
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5. Analysis of the mSUGRA parameter space 

Supersymmetry has attracted attention since it was first proposed for its many intriguing 
theoretical features and also for providing an elegant solution to the hierarchy problem 
(for an excellent introduction to the supersymmetry formalism and phenomenology see 

L3J). 

The minimal phenomenologically viable supersymmetric extension of the SM is 
called the Minimal Supersymmetric Standard Model (MSSM). The MSSM has appealing 
features: it is in much better agreement with the assumption of unification of the gauge 
couplings than the Standard Model [2T] . besides, it provides with a good dark matter 
candidate when the lightest supersymmetric particle (LSP) is the neutralino (see for 
instance [22]). In the MSSM every known particle is associated to a superpartner to 
form either a chiral or a gauge supermultiplet, whose spin components differ from each 
other by 1/2. In addition it has two SU(2) doublet complex scalar Higgs fields (H u 
and Hd) in order to give masses to the up and down type particles and to avoid gauge 
anomalies. The ratio of the vacuum expectation values of the neutral part of these two 
Higgs doublets is known as tan/5 =< H® > / < >= v u /vd- The MSSM includes a 
discrete matter parity called R parity. All the SM particles have charge +1 under this 
symmetry, whereas the supersymmetric partners have charge —1. R parity is responsible 
for the stability of the LSP. 

Since we do not observe any of the superpartners at the current laboratory energies, 
if supersymmetry exists it has to be broken. The breaking of SUSY introduces around 
~ 100 new parameters, called soft breaking SUSY terms. These are constrained by the 
limits on flavour changing neutral currents and CP-violating processes, which lead to 
the assumption of "universality" of the soft breaking terms. The universality condition 
means that the gaugino masses have a common value at the GUT scale, and the 
same assumption is made for the scalar (squark and sleptons) masses and the trilinear 
couplings respectively. After SUSY and electroweak symmetry breaking, there is mixing 
between different squarks, sleptons, and Higgses with the same electric charge, and 
between the gauginos and electroweak Higgsinos. Thus, the low energy mass spectrum 
of the MSSM contains the usual SM particles, squarks and sleptons, four neutralinos 
(denoted by Xi 2 3 4)5 two charginos (Xi 2 )j the gluino, and five Higgs bosons: the usual 
light SM-like Higgs M Higgs , a heavy neutral one M H , two charged ones M^, and a 
pseudoscalar Ma- 

Minimal supergravity (mSUGRA) is one of the better studied SUSY breaking 
scenarios [23] . In this class of models supersymmetry is broken spontaneously in a hidden 
sector that connects only through gravitational-strength interactions with the MSSM or 
visible sector. In the visible sector these interactions induce the appearance of the soft 
SUSY breaking terms. These are determined by only five parameters, a universal mass 
for the gauginos at the GUT scale m^, a universal mass for the scalars m , a universal 
trilinear coupling Aq, tan (3, and the sign of the Higgsino mass parameter fi. This drastic 
reduction in the number of parameters facilitates the scanning of interesting regions of 
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the parameter space. This is usually done by fixing two parameters, for instance Aq 
and tan (3. It is also possible to vary the four continuous parameters freely using Monte 
Carlo techniques with interesting results [Ml EH] • 

We are finally in a position to use both the abundance and entropy criteria, to 
compute the total mass density of neutralinos today, and constrain the region in the 
mSUGRA parameter space where both criteria are fulfilled. Out of the five parameters 
of the model, we will consider that the Higgsino mass parameter has a positive sign. 
This consideration is based on studies of the anomalous magnetic moment of the muon 
<7 M — 2, where SUSY models with a positive sign of \x can give a much better agreement 
with the experimental value of g^ — 2 than in the Standard Model, whereas negative 
values do not solve this problem [26j-|29j. 

Our strategy is then to explore broad regions of values for the other four parameters, 
by means of a bidimensional analysis in the mo-m^ plane with different fixed values of 
Aq and tan/3. It is important to mention that we are not presenting an exhaustive search 
in all the possible regions, but we concentrate on those regions which have received more 
attention in the literature, see for example [I]. 

In Fig. (J2j), we present the results for tan/3 = 10, and for three values of Aq, namely 
Aq = 1000, 0, —1000 GeV, shown in the top, middle and bottom panels respectively. The 
figure shows the so called bulk and coannihilation regions. The yellow region (lower right 
corner) is where the stau f is the LSP, the lighter and darker areas (red and blue for 
the online version in colors) define the allowed regions for the EC and AC respectively 
according to the observed value of Q-dm- The area of the EC region depends on the 
size of the interval of values of the parameter a, Eq. fPT]) . the lower and upper bounds 
of a determine the upper and lower boundaries of the EC region. As can be seen from 
the figure, the region where both criteria are fulfilled is very small, in fact, only for the 
highest values of a there is an intersection between both criteria. This behavior holds 
for all values of Aq in the interval [—1000, 1000] GeV; here we are showing the results 
only for the extreme and central values of aJJJ. 

Repeating the same procedure for larger values of tan/3, we find that the intersection 
region for both criteria becomes larger, getting more significant for larger values of this 
parameter. This is clearly shown in Fig. ([3]), equivalent to Fig. (j2j), but for tan/3 = 50. 
In this case the bottom panel is for Aq = —500 GeV. It is clear from the figure that 
for this value of tan/3, both criteria are consistent, there is a large intersection area for 
values of i n the interval [0, 1000] GeV. For negative values of Aq, the intersection 
region decreases as Aq does so, see the bottom panel of Fig. ([3]). For even lower values 
of Aq the intersection becomes insignificant. 

In Fig. we present the same analysis but for a different region of the parameter 
space (high values of m ), known as Focus Point region, and for the central value Aq = 0. 
The situation is consistent with the previous result, both criteria intersect for tan (3 = 50 
and there is no intersection for tan/3 = 10. 

|| For Fig. |J2) and the following Fig. (TJJ, the disconnected regions are caused by discreteness on the 
grid values chosen to explore the mp — mi m plane. 
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Figure 2. Allowed regions in the parameter space for AC (lighter gray, red) and EC 
(darker grey, blue) for the mSUGRA model with sgn^i = +, tan/3 = 10, and A — 1000 
GeV, top panel, A — GeV, middle panel, and A n = —1000 GeV, bottom panel. The 
yellow region (lower right corner) is where the f is the LSP. The figure shows the so 
called bulk and coannihilation regions. 

This analysis allows us to arrive to one of the main results of our work. The use of 
both criteria favors large values of tan/?. 

It is interesting to point out that, within the AC formulation, the dark matter 
density constraint reduces effectively the allowed regions in the mSUGRA parameter 
space to strips pointing to an almost linear trend between m Q and m^, at least in the 
coannihilation and focus point regions, as can be clearly seen in the darker (blue) areas 
in Figs. (j2l O HJ). This type of behavior has been noticed previously in [30], where the 
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Figure 3. The same as Fig. but for tan/3 = 50, and now A = —500 GeV in the 
bottom panel. 

authors even give some parameterizations for these so called "WMAP lines" . Notice that 
this type of relation is also shown within the EC description in the above mentioned 
regions, although the "WMAP lines" are not as narrow as in the AC formulation, 
because of the wide range in the a parameter Eq. (IT7|) . A deeper analysis of the 
connections among the mSUGRA parameters could reveal the original cause for such 
behavior, shown in both, quite independent, criteria. 

The regions that were compatible with the abundance and entropy criteria were 
analyzed to see what restrictions they imposed on the SUSY spectra. This was done 
just in order to show the expected mass spectra resulting from these constraints. The 
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Figure 4. Allowed regions in the parameter space for AC (lighter gray, red) and EC 
(darker grey, blue) in the mSUGRA model with Aq = 0, sgn/j = +, tan/3 = 10, top 
panel, and tan/3 = 50, bottom panel. The figure shows the so called Focus Point region. 



analysis was performed using micrOMEGAs linked to Suspect [11], using the default 
input values of micrOMEGAs, for instance: 

M top = 175GeV, (19) 
m bot = 4.23 GeV. (20) 

The Higgs bosons masses, as well as the rest of the superpartner masses, depend directly 
on the Universal gaugino mass mu2 and the Universal scalar mass mo- The first 
noticeable fact is that the bound on the Higgs mass [31] : 

mttiggs > 114 GeV, (21) 

favors the results with large tan/3 for the bulk and co-annihilation regions, as can be 
seen from Fig. which shows the Higgs mass M Higgs plotted against the mass of the 
LSP. In the case of both the bulk and co-annihilation regions, for tan f3 = 10 only a 
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small region of parameter space is allowed, with m x ~ 160 GeV. In this same case for 
tan /3 — 50 the allowed SUSY spectra starts for m x > 141 GeV. In the case of the focus 
point region the Higgs mass does not impose any further constraint, and m x > 114 
GeV. 
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Figure 5. The figure shows the Higgs mass plotted versus the LSP mass. Points 
above the dashed line are the allowed values for the LSP. 

Fig. (jSJ) top panel, shows the different Higgs bosons masses plotted against the 
pseudoscalar one Ma, for tan/3 = 50. As can be seen from the graph, the heavy neutral 
Mjj and charged Mh+ Higgs bosons are much heavier than the lightest one, and almost 
degenerate in mass. The lightest supersymmetric observable particle LOSP is plotted 
against for tan/? = 50 for the bulk and co-annihilation regions. The region above 
the dashed line shows the allowed values after the constraint on Mm ggs is taken into 
account. A comparison of the LSP with the LOSP shows that for this region, for smaller 
values of mi/2 the LOSP will be the lightest chargino X s1 (which is practically degenerate 
in mass to the second lightest neutralino), and very close in mass with the stau. As mi« 
increases the LOSP is the stau, again with values very close to the LSP, as expected. 
For the focus point region the LOSP is the lightest chargino X\i followed by the second 
neutralino x°j with almost degenerate masses. Similar results hold for tan/3 = 10 albeit 
for a much more reduced region in parameter space. The gauginos Mi, M 2 , M 3 are 
plotted against rayi in the bottom panel for tan/? = 50. The masses of these particles 
on the focus point region are very similar to the ones shown in the figure (they basically 
lie along the same lines). In both, large and small tan/? cases, the gluino is much 
heavier than the other two gauginos, as expected, and this difference increases with a 
larger value of the neutralino mass, or equivalently, with larger mi/2- 

6. Conclusions 

We have followed the novel idea of [Tj to introduce a new criterion to constrain the 
mSUGRA parameter space, which uses the assumption of entropy consistency for the 
initial and final states of a neutralino gas. Using the program micrOMEGAs, we 
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Figure 6. The lightest Muiggs and heavier Higgs bosons Mh, Mjj plotted against 
Ma for tan/3 = 50 (top panel). The lines to the right of the graph correspond to the 
focus point region, where the masses are in general heavier (except for Mniggs)- The 
mass of the LOSP plotted against rriin for tan/3 = 50 in the bulk and co-annihilation 
regions (middle panel). The masses of the gauginos Mi, M2, M3 plotted against m 1 / 2 
for tan/3 = 50 also in the bulk and co- annihilation regions (bottom panel). 

explored with precision which regions simultaneously satisfy this criterion and the usual 
abundance criterion previously used in the literature. 

Combining both, and using one of the most recent analysis that combines 



Constraining the mSUGRA parameter space using the entropy of dark matter halos 21 



observations coming from different sources to constrain the actual dark matter density, 
we were able to show that both criteria are compatible for large values of tan j3, in 
particular for tan/3 = 50, and that for lower values (tan {3 = 10), the coincidence is scarce. 
The result then favors the scenario where the ratio of the vacuum expectation values 
for the neutral supersymmetric Higgs bosons is large. Some other SUSY extensions of 
the SM, taken together with LEP data [32J, can also favor large values of tan/3 (for 
instance Grand Unified or supergravity models with Yukawa coupling unification (see 
for instance (33J, [33] ) , or Finite Unified Theories (see for instance [351 136] )). 

Moreover, we went a step further and analyze the mass spectrum of the SUSY 
particles and obtained a bound for the neutralino mass for the bulk and co-annihilation 
region, m x > 141 GeV, which is higher than the actual experimental constraint, m x > 46 
Gev [371. 



In the process to obtain the presented results, we described carefully a method 
to obtain an empirical interval on the values for a, a parameter reflecting the present 
ignorance on the appropriate statistical-mechanics description on dark matter systems. 
We used a halo model with a constant density core in the center followed by an NFW 
profile to describe dark matter dominated systems, which are the astrophysical systems 
for which our model is suitable. We found that, although a is a model dependent 
quantity, the bounds reported in this paper for its value (Eq. (fTTj) ) are reliable, 
considering the different range of observational data we used to constrain its value. 
More precise observations on the line-of-sight velocity dispersion of dwarf spheroidals 
are promising to improve our results. 

Finally, we want to remark that the analysis presented in this work is robust and it 
can be done for any other particle claimed to be candidate for dark matter, or for other 
supersymmetric extensions of the Standard Model. Also, our results can be extended for 
other regions of the parameter space, for instance, analyzing A ^ for the Focus Point 
region, fi < and exploring the "rapid annihilation funnel" region (see for example right 
panel, Fig. (4) of [38]). 
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